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The ridge in proton-proton collisions at the LHC 
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We show that the key features of the CMS result on the ridge correlation seen for high multiplicity 
events in yfs — 7 TeV proton-proton collisions at the LHC can be understood in the Color Glass 
Condensate framework of high energy QCD. The same formalism underlies the explanation of the 
ridge events seen in A+A collisions at RHIC, albeit it is likely that flow effects may enhance the 
magnitude of the signal in the latter. 



I. INTRODUCTION 
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A very recent preprint [lj from the CMS collabora- 
tion at the LHC reports the observation of a ridge-like 
structure of correlated charged particle pairs with mo- 
menta in the range p± , q± ~ 1-3 GeV in high multiplic- 
ity (with TV > 110) events in proton-proton collisions at 
y/s = 7 TeV. This ridge is a feature on the "near side" 
of the two particle correlation, around A</> w azimuthal 
separation between the two particles. It extends at least 
up to Arj R3 4.8, the limit of acceptance of the detector 
components used in the analysis. The ridge is seen only 
in the moderate p± , q±_ range and systematically vanishes 
for p±,q± < 1 GeV and p±,q± > 3 GeV. A plot summa- 
rizing this structure in the CMS data is shown in fig. [TJ 



A similar ridge-like structure was previously seen in 
nucleus-nucleus collisions at RHIC @t3|- The ridge was 
seen in high multiplicity (central) events in Cu+Cu col- 
lisions at </s — 62.4 GeV and in Au+Au collisions at 
y/s = 200 GeV. The STAR detector observed this cor- 
relation for both p^-triggered @, |j| and untriggered Q 
pair correlations in the whole STAR TPC acceptance of 
A77 < 2. The PHOBOS experiment [6j observed the p±- 
triggered correlation at much larger rapidity separations, 
initially up to A 77 ~ 4, extended more recently Q to 
A77 ~ 5. The long range correlation structure disap- 
pears for lower multiplicity peripheral events in nucleus- 
nucleus collisions and are also absent in deuteron-gold 
and proton-proton "control" experiments at the same en- 
ergies. The STAR and PHOBOS correlation results are 
shown in fig. rj] 



where Tf. is the freeze-out time of the particles. This 
causality argument is illustrated in fig. [3J If we assume 
this freeze-out time to be of the order of a few fermis, a 
rapidity correlation of A77 ~ 5 units suggests that these 
correlations must be formed nearly instantaneously after 
the collision or must preexist in the incoming projectiles. 
They are therefore very sensitive to the strong color fields 
present in the initial stage and their correlations. 

In a recent paper (see also [13, El), it was shown 
in detail how these long range correlations arise as a 
consequence of the saturation of gluons with momenta 
k± < Q s in the nuclear wave functions before the col- 
lision. Here Q s > Aqcd is a semi-hard saturation 
scale [HI, EH , which provides a measure of the range of 
color correlations in a nucleon or nucleus. This scale 
is large for either high energies or large nuclei or both, 
suggesting the presence of strong universal dynamical 
color correlations at distances smaller than the scale of 
confinement dynamics 1/Aqcd in QCD. Because Q s is 
large, a quantitative understanding of the highly non- 
perturbative dynamics of saturated gluons in feasible in 
weak coupling, and is realized in the Color Glass Conden- 
sate (CGC) effective theory [H-Ei of hadron and nuclear 
wavefunctions at high energies. The reference [3] pro- 
vides a quantitative comparison with the PHOBOS data 
and makes predictions for similar correlations in A+A 
collisions at the LHC. 

When two hadrons or nuclei collide at high energies, 
the CGC framework predicts that strong longitudinal 
chromo-electric and chromo-magnetic fields arc produced 



Why is the ridge interesting? One can show from a 
simple argument based on causality Q that if long range 
rapidity correlations between particles exist, the correla- 



This argument assumes that a produced particle's momentum 
space rapidity is tightly correlated with its space-time rapidity. 
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FIG. 1: Top: 3-D display of the two particle correlation R 
as a function of A77 and A<f> for minimum bias and high mul- 
tiplicity events in two different p± windows. Bottom: The 
associated yield in different p± windows as a function of the 
number of charged particle tracks. From ref. [J. 

that are nearly boost-invariant (l7l [Tp ; this form of mat- 
ter has been called the Glasma We should em- 
phasize here at the outset that the properties of the 
Glasma are not computed in an ad hoc model but fol- 
low from high energy factorization theorems which re- 
late multi-particle dynamics in the Glasma to multi- 
parton correlations in saturated hadron/nuclear wave- 
functions [ll|, HO, [2l[ . A key feature of the Glasma is that 
the transverse correlation length of these strong longitu- 
dinal chromo-electric and chromo-magnetic fields is 1 /Q s . 
Two-particle [1, [22|, HH and three-particle correla- 
tions were computed in this Glasma "flux tube" picture. 
We should note that the n-particle correlations can be 
shown to satisfy a negative binomial distribution [2o|- 
providing a microscopic derivation of this widely used 
multiplicity distribution. These correlations are inde- 
pendent of Arj up to quantum corrections that become 
important for Ar/ > l/a s . In heavy ion collisions the 
collimation in A(j> comes primarily from the subsequent 
radial flow of the correlated particles (26l - [3l| . 

However, as discussed in more detail below, an intrinsic 



FIG. 2: Top: 3-D display of the two particle correlation 
(analogous to R in fig. [TJ for px-triggered events from the 
STAR collaboration |5| (note that the away-side peak around 
Acj> = n seen in fig. [T] has been removed). Bottom: Two- 
particle correlation data from the PHOBOS collaboration [f| 
that reveals a long range component. The curves shown are 
obtained by adding our result (in eq. @) to the short range 
correlation from PYTHIA. 



A(j) collimation is also present independently of the effects 
from flow in the later stages of the collision. It is very 
weak for the small Q s values that are attained in p+p, 
d+A and peripheral A+A collisions at RHIC. At LHC 
energies, for the large multiplicity cuts performed by the 
CMS experiment (which select central impact parame- 
ters), the signal is stronger because Q s is larger but not 
as strong as in central A + A collisions, because the effect 
of radial flow is smaller or absent. Indeed, "predictions" 
previous to the CMS announcement existed [321 ] but were 
not submitted for publication because it appeared in- 
conceivable that such a small signal would be detected 
experimentally-this places the tour de force measurement 
by CMS in perspective! 

We shall now show some of our results that exhibit 
the same qualitative features as the CMS results. More 
quantitative comparisons would require a more detailed 
understanding of the impact parameter dependence af- 
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FIG. 3: This figure illustrates the argument from causality 
that long range correlations of particles (denoted A and B 
here) must occur at very early proper times. The doubly 
shaded region, corresponding to the intersection of the light 
cones of the two particles, is the space-time location where 
correlations can be formed. 



fccting the normalization, experimental efficiency and ac- 
ceptance effects and the contribution of radial flow. We 
hope to be able to make more detailed comparisons in 
the near future. 
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FIG. 4: The unintegrated gluon distribution as a function of 
transverse momentum squared for running coupling for three 
values of the momentum fraction x. 



II. LONG RANGE TWO-PARTICLE 
CORRELATION 



The strength of two-particle correlations is conve- 
niently represented by the function 
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The expression in eq. @ corresponds (qualitatively) 
to the quantity R(A<fi, Arj) plotted in the CMS paper. 
Our result for the correlated two-particle distribution in 



eq. ([2]) can be expressed as [9[ 
dN 2 a 2 
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d 2 p ± dy p d 2 q ± dy ? 16^ 10 (iV 2 - l) 3 p 2 ± q 2 ± 

x J d 2 k ± {<$> 2 A {y p ,k ± )<S> B (y p ,p ± -k x ) 

x [$ s (t/ 9 ,qx +k x ) + $B(j/«,q_L - k x )] 
+$|(y 9 ,k ± )$ A (y p ,p ± - kx) 

x [$ A (y g ,q ± +kx) + &A(y g , q± - kx)]} • 

Here the ^'s are unintegrated gluon distributions (UGD) 
per unit of transverse area in projectiles A and B and 
S± is the transverse overlap area of the two hadrons. 
Eq uation [3] is based on the formalism developed in [ToL 
[TT| and is derived in ref. @. The single inclusive gluon 
spectrum can be expressed in the same notation as (33" 
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The above expressions (eqs. [3]and[4]) are valid to leading 
logarithmic accuracy in x and for momenta p±, q± > Qs- 
The important ingredient in the above expressions 
is the unintegrated gluon density which is a universal 
quantity and can be constrained from fits to DIS on 
hadron [3(| and nuclear targets @-these unintegrated 
gluon densities have been used to compute single inclu- 
sive and double inclusive distributions in d+Au collisions 
at RHIC (3?| • The only difference between different tar- 
gets are in the initial conditions determined at an initial 
xo-chosen to be xo = 0.01. The evolution of the UGD 
with rapidity is controlled by the Balitsky-Kovchegov 
(BK) equation (38l . [39j . which is a non-linear evolution 
equation describing both gluon emission and multiple 
scattering effects. We note that our unintegrated gluon 
distributions have recently been shown to be equivalent 
to those derived in the framework of transverse momen- 
tum dependent (TMD) parton distributions [40j. In fig. 0] 
we show the structure of this unintegrated gluon dis- 
tribution and its evolution with rapidity for a running 
coupling. The running of the coupling slows the evolu- 
tion significantly which is essential to describe DIS and 
hadron scattering data. 



III. THE RIDGE IN PROTON PROTON 
COLLISIONS 

In proton-proton collisions, unlike nucleus-nucleus col- 
lisions, it is not clear that there can be a sufficient amount 
of transverse flow to provide the collimation around 
Acj) w necessary to explain the ridge. Blast wave fits to 
proton single particle spectra at low p± are suggestive of 
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non-vanishing radial flow, and this in principle could pro- 
vide an additional collimation of the signal. We also note 
that in the highest multiplicity proton-proton events the 
charged particle multiplicity per unit rapidity is compa- 
rable to that found in semi-central Cu-Cu collisions. The 
possibility that additional collimation of the signal due 
to flow might be needed to reproduce the full strength 
of the measured correlation cannot be completely ruled 
out. However, the absence of a correlation signal for small 
transverse momenta in proton-proton collisions (contrary 
to the nucleus- nucleus case) suggests that hydrodynami- 
cal flow would not be the dominant contribution here. 



* 8 



3 



8050? 



3000550000m 

P, Up \ 



<h Vq 

>ooooooooooo@ 550C5 



P -k 



FIG. 5: A typical diagram which gives an angular collimation. 

The main point of this paper is that there is an in- 
trinsic angular correlation (in addition to the long range 
rapidity correlation) coming from the particle production 
process on transverse distance scales 1/Q S much smaller 
than the proton size. One of the diagrams that gives 
a collimation for A(f> w is shown in fig. [5] There is 
only a single loop momentum kj^ in these two-particle 
production diagrams, as opposed to uncorrelated produc- 
tion. The transverse momenta flowing into the "blobs" 
for hadron A in this diagram are |pj_ — kj_ | and | qj_ — kj_ | , 
respectively. If we combine this with the fact that the 
unintegrated distributions peak about Q s (fig. 0} it fol- 
lows that the largest contribution is obtained when both 
| — kj_ | ~ Q s and |q_i_ — kjj ~ Q s . These conditions are 
satisfied simultaneously, for one and the same kj_, if pj_ 
and are parallel, leading to angular collimation. The 
scale of the angular dependence is 0(1). Also, if para- 
metrically p± and q± are much smaller or much greater 
than Q s , the collimation disappears. 

In contrast, in collinear factorization, at leading order 
both gluons are produced from the same ladder leading to 
collimation about A<fi ~ ir ("dijet"). One may therefore 
expect that diagrams such as fig. O give a large contribu- 
tion to C(p_L,t/ p ,q_L,y g ) when Acj) < 7r and \y p — y q \ > 1 
(see, also, discussion in ref. fioj)- 

Our results are summarized in figs. [6] and [7] In fig. [6] 
we show C2(p_l, y p , qxj Vq) for p± s» q± as a function of 
Qq, the saturation scale at the initial condition of the 
evolution x — Xq. A larger saturation scale corresponds 
to a more central collision and thus to a higher multiplic- 
ity. Thus our result shown in fig. [5] agrees with the CMS 



observation of the ridge correlation becoming visible for 
higher multiplicity events. The exact relation between 
the value of Q s and the observed multiplicity depends 
on the overlap area of the two protons and on the im- 
pact parameter dependence of the saturation scale (see 
e.g. [ill, [42j]). Establishing the normalization without ad- 
ditional free parameters would require going beyond the 
&X-factorized expression (jj) for the single inclusive spec- 
trum (43|. However, in the appropriately normalized cor- 
relation ((2J these factors cancel and our result does not 
depend strongly on the size of the overlap area. The re- 
sult shown in fig. [5] is obtained by integrating over A<j) 
for < A<j) < 7r/2. For a fixed p±, we see that the 
correlation increases with centrality. The systematics is 
striking. For all centralities, the collimation is weak for 
low p±. For p± > 1 GeV, the correlation grows but at 
p± > 6 GeV, the correlation decreases below the value 
at p± — 0.5 GeV and is smaller for all higher p± values. 
The reader should keep in mind that this result is for 
gluons. Fragmentation effects may approximately cancel 
in the C2 ratio but the p± of the hadrons is of course 
lower than that of the gluons. While this plot cannot 
be compared directly to fig. [T] (bottom) it shows a non 
trivial systematics very similar to it. 

In fig. [7] (top) we plot the p± dependence of C2 . We 
see that the distribution is peaked approximately at the 
value of Q s evolved to the rapidities of interest and that 
the correlations drop off sharply at large p± . For low p± , 
one has less theoretical control on how the distributions 
fall off but the trend is unmistakable. In fig. [7] (bottom), 
we show the rapidity dependence for fixed p± = q± = 
2 GeV as a function of y q , given two different y p and 
two different values of Qq. For y p = 0, we see that the 
distributions are symmetric around y q = 0. This is not 
the case for y p = —3, because QCD evolution (gluon 
radiation) between the projectiles to the produced gluons 
is asymmetric. The systematics could in principle be 
investigated to test our predictions in detail. 



IV. SUMMARY 

The physics of gluon saturation opens up a novel do- 
main of many body QCD at high energies and may con- 
tain many surprises which provide deeper insight into 
the fundamental properties of strong interactions. One 
of these is the ridge in proton-proton collisions and it is 
likely that many more will be revealed at the LHC. The 
properties of saturated gluons, in particular the evolu- 
tion of multi-parton correlations with energy, can be com- 
puted in the Color Glass Condensate effective field the- 
ory. High energy factorization theorems allow us to relate 
these multi-parton correlations in the hadron and nuclear 
wavefunctions to inclusive multi-parton final states in col- 
lisions. The key expression eq. ([3]) is one of the results 
of this approach. It provided a detailed understanding of 
the ridge in A+A collisions at RHIC. It predicts a ridge 
in p+p collisions. 
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FIG. 6: The correlation function C2 for different p± as a function of Qq, the saturation scale at the initial condition of the 
evolution x = xq. Increasing values of Qo correspond to smaller impact parameters and larger total charged multiplicities. 



We computed in this paper the two-particle correla- 
tion C2 (eq. ©) and showed that it shares the same key 
features as the quantity R measured by the CMS exper- 
iment: 

• It is long range in rapidity, and exhibits a collima- 
tion for high multiplicity events for a narrow win- 
dow of a few GeV in the transverse momenta of 
the pair particles. The effect becomes systemati- 
cally weaker both below and above the kinematic 
window as seen in the CMS result. 

• The two-particle correlation has the same strength 
for both like and unlike sign pairs. This feature was 
seen already at RHIC and is consistent with gluon 
emission from a source that is (nearly) uniform in 
rapidity "the glasma flux tube" . It is not consistent 
with emission from jets. 

• The two-particle correlation is relatively flat in the 
771 versus 772 plane. Again, this feature is natural for 
gluon emission from nearly boost invariant classical 
sources. We predict a similarly flat signal strength 
in three-particle correlations [24[ . 

Correlations of trigger particles at non central rapidities 
with associated particles at other rapidities can provide 
sensitive tests of this picture. In order to make a quanti- 
tive comparison one needs to add the short range jet com- 
ponent (which can be generated with one of the available 



event-generators) to our long range correlation. After 
including fragmentation one could then form the same 
quantities as measured in experiment. This is left for 
further study. 
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